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ABSTRACT 

We compute the cross correlation of the intensity and polarisation from the 5-year WMAP 
data in different sky-regions with respect to template maps for synchrotron, dust, and free- 
free emission. We derive the frequency dependence and polarisation fraction for all three 
components in 48 different sky regions of HEALPix (iVgidc — 2) pixelisation. The anomalous 
emission associated with dust is clearly detected in intensity over the entire sky at the K 
(23 GHz) and Ka (33 GHz) WMAP bands, and is found to be the dominant foreground at low 
Galactic latitude, between b = —40° and b = +10°. The synchrotron spectral index obtained 
from the K and Ka WMAP bands from an all-sky analysis is f3 — —3.32 ± 0.12 for intensity 
and /3 = —3.01 ± 0.03 for the polarised intensity. 

The polarisation fraction of the synchrotron is constant in frequency and increases with 
latitude from « 5% near the Galactic plane up to « 40% in some regions at high latitude; 
the average value for |6| < 20° is 8.6 ± 1.7(stat) ± 0.5(sys)% while for \b\ > 20° it is 
19.3 ± 0.8 (stat) ± 0.5 (sys) %. Anomalous dust and free-free emission appear to be relatively 
unpolarised. Monte carlo simulations showed that there were biases of the method due to 
cross-talk between the components, at up to « 5 % in any given pixel, and « 1 . 5 % on average, 
when the true polarisation fraction is low (a few per cent or less). Nevertheless, the average 
polarisation fraction of dust-correlated emission at K-band is 3.2 ± 0.9 (stat) ± 1.5 (sys) % or 
less than 5 % at 95 % confidence. When comparing real data with simulations, 8 regions show 
a detected polarisation above the 99th percentile of the distribution from simulations with no 
input foreground polarisation, 6 of which are detected at above 2a and display polarisation 
fractions between 2.6 % and 7.2 %, except for one anomalous region, which has 32 ± 12 %. 
The dust polarisation values are consistent with the expectation from spinning-dust emission, 
but polarised dust emission from magnetic-dipole radiation cannot be ruled out. Free-free 
emission was found to be unpolarised with an upper limit of 3.4 % at 95 % confidence. 

Key words: ISM: general - Galaxy: general - cosmology: diffuse radiation - cosmology: 
cosmic microwave background - radio continuum: ISM 



1 INTRODUCTION 

During the last few years, great advances have been made in mea- 
suring the Cosmic Microwave Background (CMB) anisotropies. 
More advances are expected from the Planck mission that will pro- 
vide all-sky temperature and polarisation CMB measurements with 
high accuracy and unprecedented frequency coverage. As CMB ex- 
periments reach higher sensitivity, it becomes necessary to better 
understand and remove the contribution from competing Galac- 
tic and extragalactic astrophysical signals in the same frequency 
bands. It is well known that the emission of these foregrounds ri- 
vals the temperature signal over a significant fraction of the sky 
and dominates the polarisation signal over most of the sky (Kogut 



et al. 2007). As the temperature power spectrum has already been 
measured with great precision (e.g., WMAP5, Nolta et al. 2009; 
ACBAR, Reichardt et al. 2009; Boomerang, Jones et al. 2006; CBI, 
Sievers et al. 2009; VSA, Dickinson et al. 2004; QUaD, Pryke et 
al. 2009), most observations are now concentrating on polarisation 
measurements (e.g., BICEP, Chiang et al. 2009; QUaD, Brown et 
al. 2009). For the study of polarisation anisotropies however, fore- 
grounds may be the limiting factor, as the primordial signal is much 
weaker than intensity while the foregrounds are not well charac- 
terised at the relevant frequencies at this time. For this reason, suit- 
able component separation techniques have been developed and 
are being applied to current data (Kogut et al. 2007; Eriksen et 
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al. 2008; Gold et al. 2009; Dunkley et al. 2009; Delabrouille & 
Cardoso 2009), in preparation for Planck data, and possible future 
missions (Leach et al. 2008; Betoule et al. 2009). In some instances 
such techniques rely on prior information acquired on the emission 
of a specific foreground, for which a prior characterisation is neces- 
sary. Moreover, the study of foregrounds in the microwave band is 
of astrophysical interest on its own, as it can be informative about 
physical phenomena occurring in our Galaxy and in others. 

Foregrounds are either Galactic or extragalactic in nature. In 
this paper we are interested in the Galactic emission, which dom- 
inates over extragalactic emission at intermediate to large scales 
QL 1°)- The Galactic emission has three main components: dust, 
synchrotron, and free-free, of which only dust and synchrotron are 
expected to show significant polarisation. In contrast with extra- 
galactic emission, we do not expect the Galactic foreground to be 
isotropic. For this reason. Galactic foregrounds need a characteri- 
sation that is both scale dependent and position dependent. 

The most common methods for extracting foreground infor- 
mation from CMB data involve fitting a given theoretical model to 
the data. These fitting procedures typically assume that a specific 
foreground follows a power-law in frequency and is characterised 
by a given polarisation fraction and a (fixed) angle. For instance, 
this is the approach of the WMAP teams (Kogut et al. (2007) for the 
WMAP 3-year and Gold et al. (2009) for the WMAP 5-year data). 

By so doing, however, the diverse behaviour of different fore- 
grounds in different regions of the sky may be neglected. For in- 
stance, a spinning-dust contribution may not have a power-law fre- 
quency dependence, and may make different contributions in dif- 
ferent areas of the sky. A variety of dust grains may populate dif- 
ferent areas of our galaxy, and have different contributions to the 
intensity and polarisation signal at the observed frequencies. The 
synchrotron emission clearly depends on the energy spectrum of 
electrons in a given region, while the polarisation signal may be af- 
fected by depolarisation effects which are more severe in denser or 
less uniform regions. 

In this paper, we aim to characterise the contribution of the 
three Galactic components in the microwave band in temperature 
and polarisation by studying the WMAP5 data (Nolta et al. 2009; 
Hinshaw et al. 2009), including their frequency and spatial depen- 
dencies. To achieve this goal we develop a foreground analysis that 
assumes no frequency model and only relies on a morphological 
characterisation of foregrounds as derived at non-CMB frequen- 
cies. Specifically, we use a cross-correlation (C-C) analysis with 
standard foreground templates. Such an analysis method has been 
widely used by other authors studying CMB intensity data (e.g., 
de Olivera-Costa et al. 1999, 2002; Banday et al. 2003; Bennett et 
al. 2003; Davies et al. 2006), but has not been applied to polarisa- 
tion data. An advantage of using such a method is the possibility of 
detecting particular emission processes associated with one of the 
three foreground templates even if its frequency and/or polarisation 
characteristics are unknown. For example, this may be true in the 
case of a spinning-dust component. 

A possible limitation of the C-C method is that it allows char- 
acterisation of only foreground components whose spatial distribu- 
tion is well traced by the templates. These templates were derived 
at very different frequencies than the CMB measurements, and may 
not fully represent the morphology of a particular component in the 
range of frequencies at hand. Most analyses of WMAP5 polarisa- 
tion data are based on modelling the various emission components 
in the Q and U maps; here we will perform our analysis in the po- 
larisation intensity P. This choice avoids the problem of modelling 
the Galactic magnetic field to compute the polarisation angle, or 



assuming a specific polarisation angle as a function of position for 
each component. Under the hypothesis that all components have 
the same polarisation angle (even if possibly a different polarisa- 
tion intensity), the method used in analysing the temperature maps 
can be readily extended to the polarisation maps. Although a more 
general approach may be sought, as both dust grains' orientation 
and synchrotron emission are dictated by the direction of the mag- 
netic field, this assumption may not be unreasonable. 

Section|2]describes the WMAP data used in this work and the 
choices made about spatial resolution and masking, while Section[3] 
gives an introduction to the various sources of foreground emission 
and the spatial templates that will be used to trace each component. 
The cross-correlation analysis technique and the computation of 
uncertainties are derived in Section|4l The results of the analysis in 
both temperature and polarisation are presented in Section[5] along 
with separate discussions for each of the 3 emission components. 
The main conclusions are summarised in Section[6l 



2 WMAP DATA 

We use the WMAP 5-year data (Hinshaw et al. 2009), consisting 
of 5 full-sky maps at frequencies of 22.8 GHz (K-band), 33.0 GHz 
(Ka-band), 40.7 GHz (Q-band), 60.8 GHz (V-band), and 93.5 GHz 
(W-band), as provided by the LAMBDA website^ . These maps are 
provided in HEALPix^ format at a resolution of Neidc = 512 
(number of pixels on the sphere = 12 x N^nc)- The maps are then 
downgraded to a HEALPix resolution of A'aide = 32, using the 
ud.grade'^ HEALPix routine, to give a total of 12,288 pixels. We 
chose this resolution as a compromise between obtaining a rea- 
sonable signal-to-noise ratio and having small enough pixels to 
keep the information about foreground variations needed in the C- 
C analysis; with A'^sidc = 16 the resolution is becoming too coarse. 
Larger pixels are also likely to invalidate the assumption of a con- 
stant polarisation fraction. Parallel transport of polarisation vectors 
on the sphere is not considered in the averaging since the effect is 
at the 10^* level or lower. 

We convert the data from thermodynamic temperature to 
brightness (antenna) temperature using the conversion: 



where x = {hu)/{kTQ) and To is the CMB temperature of 2? 725 K 
(Mather et al. 1999); this corresponds to a correction of 1% at K- 
band and 25% at W-band. 

We apply the WMAP team mask KQ85 (Gold et al. 2009) to 
the WMAP maps, in order to avoid the Galactic plane and bright 
point sources. The mask is provided in a A'^sidc ~ 512 HEALPix 
resolution, so we downgrade it to our working resolution (A'^side = 
32) and consider in the analysis only the coarse pixels that are com- 
pletely outside of the A'^sido = 512 masked region. This procedure 
reduces to 8099 the number of data pixels (66% of the sky) used in 
the analysis. 

In the polarisation analysis we compute the total polarisation 
as P = ^/Q'^ + U"^. The mask and the resolution chosen are the 
same as the intensity analysis. This allows direct comparison of the 
C-C coefficients fit for the temperature and polarisation data and 
calculation of the fractional polarisation at the same sky positions. 

^ http://lambda.gsfc.nasa.gov/ 

^ http://www.eso.org/science/healpix/ (Gorski et al. 2005) 

Simple averaging is used, rather than a noise-weighted average. 
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3 GALACTIC COMPONENTS EMISSION AND 
TEMPLATES 

3.1 Synchrotron 

Synchrotron emission arises from accelerating cosmic ray electrons 
spiralling in the Galactic magnetic field. The emission depends on 
the energy spectrum of the electrons and on the intensity of the 
magnetic field, which results in significant spatial variations in the 
spectral index on the sky. The ssmchrotron spectrum is approxi- 
mated by a power law T ~ u'^, typically with /? ~ —2.7 at radio 
wavelengths, steepening to P ~ —3.0 at WMAP frequencies and 
with typical spatial variations of ±0.2. Moreover, it is known that 
a pure synchrotron spectrum steepens with frequency due to the ef- 
fects of spectral ageing, although in reality, a flattening of the spec- 
trum can occur due to the presence of multiple components (e.g., 
Kogut et al. 2007). 

Most of the information we have for synchrotron emission 
comes from low frequency radio surveys, where synchrotron dom- 
inates the sky. In particular, large-area radio surveys at 408 MHz 
(Haslam et al. 1981, 1982), 1420 MHz (Reich & Reich 1986) 
and 2.3 GHz (Jonas et al. 1998) give us a clear picture of the 
synchrotron component and these have been used extensively for 
foreground subtraction and previous cross-correlation studies. The 
drawbacks of using such templates are i) there is a large extrap- 
olation in frequency from ~ 1 GHz to WMAP frequencies that is 
likely to result in some distortion of the morphology, ii) baseline 
(offset) issues result in striping in the maps, and iii) discrete extra- 
galactic sources contaminate the diffuse emission. Several efforts 
have been made to reduce the effects of striping and source con- 
tamination (e.g., Davies et al. 1996; Platania et al. 2003) including 
the NCSA version of the Haslam et al. 408 MHz map, available at 
the LAMBDA website. We use this as our template for synchrotron 
emission. 

Synchrotron emission is, by nature, highly polarised. For a 
power-law distribution of electron energies N{E) oc E~''' prop- 
agating in a uniform magnetic field, the resulting emission is po- 
larised with fractional linear polarisation f = {p + l)/(p + 7/3) 
aligned perpendicular to the magnetic field (Rybicki & Lightman 
1979). The frequency dependence of synchrotron emission is also 
related to the electron energy distribution, T(i^) oc with spec- 
tral index /3 = — (p -|- 3) /2. For x —3,the sjmchrotron emission 
can have a maximum fractional polarisation as high as fs = 0.75. 
However, line-of-sight and beam averaging effects will tend to re- 
duce this, with typical values at high latitude of ~ 10^0%. 



3.2 Dust 

The Galactic foreground at far-infrared and submillimetre wave- 
lengths (f > 100 GHz) is dominated by thermal emission from 
warm (T ~ 10 - 100 K) interstellar dust grains. The intensity 
spectrum peaks in the range 100 - 200 [im and is well modelled 
by an emissivity-modified gray-body of the form i/^ Bv[T), where 
(T) is the Planck function at frequency v and temperature T. 
Using data at 100 and 240 i^m, Finkbeiner, Davis & Schlegel 
(1999) modelled the all-sky emission with two components hav- 
ing mean temperatures of 9.4 K and 16 K and spectral indices of 
13 = 1.7 and 2.7, respectively. The predictions at longer wave- 
lengths are in good agreement with observations at 353 GHz by 
Archeops (Ponthieu et al. 2005) and at 94 GHz by WMAP (Bennett 
et al. 2003). We use the predictions at 94 GHz, which are available 
on the LAMBDA website, as our template for dust emission. 



Polarisation from thermally emitting dust is due to aspherical 
grains whose spin-axes have become aligned locally with interstel- 
lar magnetic fields (e.g., see reviews by Lazarian 2003, 2007). The 
most efficient direction for emission and absorption/extinction is 
the long grain axis (perpendicular to the spin axis). The observed 
polarisation direction is then perpendicular to the aligning mag- 
netic field in the case of grain emission, but parallel to the field 
in the case of background-starlight extinction. The expected size of 
the polarisation depends on a mmiber of unknown factors including 
the efficiency of the alignment mechanism and the inclination of 
the magnetic field to the line-of-sight. Typical background- starlight 
polarisation has values in the ^ 1^ % range (e.g., Heiles 2000; 
Fosalba et al. 2002), while polarised emission in dense Galac- 
tic clouds have been observed at 0.5 - 10 % (e.g., Dotson et al. 
2000, 2010; Matthews et al. 2009). WMAP polarisation observa- 
tions at 94 GHz (Page et al. 2007; Kogut et al. 2007) and Archeops 
at 353 GHz (Ponthieu et al. 2005) measure dust polarisations of 
~ 1% in the Galactic plane, increasing to a few percent at higher 
latitudes. Studies of the frequency dependence of polarised emis- 
sion have found that the spectrum has a polarisation minimum at 
A ~ 350 |J.m, increasing to longer wavelengths (e.g., Hildebrand 
et al. 1999). However, these studies are limited to A = 60 - 1300 
|i.m and to bright/dense Galactic clouds. Extrapolating this spec- 
trum to the conditions expected to prevail in the high-latitude dif- 
fuse ISM, one expects a more featureless spectrum for wavelengths 
a few times longer than about 1 mm (Hildebrand & Kirby 2004). 

Other possible sources of polarised dust emission at mi- 
crowave frequencies are electric-dipole emission from "spinnng 
dust" grains, as well as magnetic-dipole emission from vibrat- 
ing "magnetic dust" grains (e.g., Lazarian & Finkbeiner 2003; 
Draine & Lazarian 1998a,b,1999; Lazarian & Draine 2000; AH- 
Haimoud, Y, Huata, C. M., & Dickinson, C. 2009). Evidence for 
the existence of spinning-dust exists in the strong correlations ob- 
served between the spatial distribution of thermal dust emission 
and that of anomalous emission at frequencies of 20 - 60 GHz. 
This anomalous emission cannot be accounted for by models us- 
ing only the standard components of free-free, synchrotron, and/or 
thermal dust emission. This dust-correlated anomalous emission is 
observed on both large-scales (e.g., Kogut et al. 1996; Finkbeiner 
2004; Finkbeiner, Langston, & Minter 2004; Davies et al. 2006; 
Hildebrandt et al. 2007; Miville-Deschenes et al. 2008) as well 
as in pointed observations of specific Galactic dust clouds (e.g., 
Finkbeiner 2004; Watson et al. 2005; Casassus et al. 2006; Dick- 
inson et al. 2009a, 2010; Scaife et al. 2009; Planck collaboration 
2011). 

The polarisation amplitude of spinning-dust emission is likely 
to be small. Lazarian & Draine (2000), modelling the alignment of 
very small dust grains in the ISM, predict that spinning dust will 
be polarised by no more than ~7% at 2 GHz, falling to < 0.5% 
above frequencies of 30 GHz. Empirically, UV observations of 
background-starlight polarisation (caused by the same very small 
grains) indicate that these grains are not as well aligned with 
the local magnetic field as are larger grains (e.g., Martin 2007). 
Magnetic-dust is expected to be polarised at much higher levels, as 
high as 40% (Draine & Lazarian 1999). The only polarisation ob- 
servations of the anomalous emission to date have found low levels. 
Battistelli et al. (2006) report 3.4li[ y% towards the Perseus molec- 
ular cloud at 1 1 GHz while Mason ct al. (2009) report an upper limit 
of 3.5% (99.7% confidence level) towards the dark cloud Lynds 
1622 at 9 GHz. Lopez-Caraballo et al. (201 1) used WMAP data to 
constrain the polarization towards the Perseus molecular cloud to 
be less than 1.0, 1.8 and 2.7 % at 23, 33 and 41 GHz, respectively. 



© 0000 RAS, MNRAS 000, 000-000 



4 Macellari et al. 



These observations are consistent with polarisation from spinning- 
dust, but cannot rule out the existence of vibrating magnetic-dust. 

3.3 Free-free 

Free-free (bremsstrahlung) emission is due to electron-electron 
scattering from warm (T^ ~ 10'* K) ionised gas in the interstel- 
lar medium. In the optically thin regime, and at radio frequencies, 
the brightness temperature is given by 



n = 8.235 X W~aT; 



=j.gh'(1 + 0.08)(EM), 



where Te is the electron temperature, EM is the emission measure, 
and a is a factor close to unity. The electron temperature in the 
warm ionised medium ranges from ~ 3000 to ~ 20, 000 K but is 
~ 8000 K at the Solar galacocentric distance. The spectrum follows 
a well defined power-law with a spectral index of —2.1 with only a 
small variation in temperature. At WMAP frequencies, the effective 
free-free spectral index is ~ —2.14 for Te ~ 8000 K. 

The well-determined spectrum can be used to separate free- 
free from other diffuse radio components, most notably, syn- 
chrotron (which is almost always steeper). However, due to the 
complexity of the components, and limited frequency range of 
WMAP, a cleaner separation can be achieved by using an appropri- 
ate spatial template. The free-free brightness temperature is propor- 
tional to EM = J rildl, i.e., the integrated colunm density squared 
along the line of sight. Fortunately, this is also true of recombina- 
tion lines such as the bright optical line. Ha. Maps of Hq there- 
fore give an almost one-to-one tracer for radio free-free emission 
(Dickinson et al. 2003) and they are now available for the full-sky 
(Finkbeiner 2003). There is some dependancy on Te and the local 
conditions, but these are relatively small. The main difficulty arises 
at low Galactic latitudes (typically |6| < 5°) where dust absorp- 
tion can reduce the Ha intensity by a large factor, thus rendering it 
useless for predicting radio free-free brightness. However, outside 
the Kp2 mask, almost all lines of sight have a small to negligi- 
ble maximum correction factor; typical values are <^ 1 mag. The 
correction factor can be estimated from the IRAS 100 /im reden- 
ning maps (Schlegel et al. 1998). We therefore use the Finkbeiner 
(2003) composite Ha map as our standard free-free template, but 
also compare a similar map produced by Dickinson et al. (2003). 

Free-free emission is intrinsically unpolarised because the 
scattering directions are random. However, a secondary polarisa- 
tion signature can occur at the edges of bright free-free features 
(i.e., H II regions) from Thomson scattering (Rybicki & Lightman 
1979). This could cause significant polarisation (~ 10%) in the 
Galactic plane, particularly when observing at high angular resolu- 
tion. However, at high Galactic latitudes, and with a relatively low 
resolution, we expect the residual polarisation to be < 1%. 



where Itot(i^, n) is the total emission at frequency u in the direc- 
tion (or map pixel) h, Ti{voi^n) is the emission from the i*'' tem- 
plate at the template frequency foi (hereafter simply Ti{h)), and 
6i{v,n) is the spectral shape of the emission template, normalised 
in every pixel such that Oili/oi) = 1. The C-C technique is ap- 
plied to sky-regions containing different pixels and assumes that 
the proportionality coefficients between the total intensity and each 
foreground component are constant within each region. That is, the 
6i{u,n) are independent of position within each region, although 
they can vary from region-to-region: 



I 



(4) 



For the emission mechanisms studied in this work we have 
n = 3 templates, although the C-C method can generally be ex- 
tended to any number n. For any given WMAP frequency u, an es- 
timate of the parameters 6i{i') can be computed by minimising x^, 
as explained in Davies et al. (2006). In the component formalism 
of Equation Q we have 



(5) 



where the sum over n is the sum over all pixels in sky-region R, Aij 
is the 3 X 3 template matrix [Ti{n)'^ ■ M^^{iy,h) ■ Tj{h)], 
and Msn(j', n) is the covariance matrix of the expected noise of 
the WMAP data and CMB signal (described in the next section). 
The objective is to invert the covariance matrix and compute the 
coefficients 8i{v) for each of the 3 Galactic foreground compo- 
nents. 

Within an area R of the sky, the C-C analysis relies on the mor- 
phological characteristics of a given component. The mean value 
of the templates' emission within R (which can be thought of as a 
constant offset in the template's map) does not provide information 
for the purpose of recovering the signal relative to each individual 
component, as it has no morphological structure. Moreover, keep- 
ing these mean offsets in the templates fixed will bias the results of 
the C-C. Since all templates, as well as all five WMAP maps, may 
have such uniform offsets they can each be represented as a single 
uniform term in the sum of Equation l[3). That is, each template and 
dataset have the form 



T{n) = T^n) ~ S[, 



(6) 



where T' is an existing template map or WMAP data, and T is its 
true value in which all constant-offsets have been removed. Equa- 
tion ^ is then re-written as 



e,{iy)T!{n)+5tot{i^) 



(7) 



where 



4 THE CROSS-CORRELATION ANALYSIS 
4.1 The Method 

The C-C method (described in Davies et al. 2006) assumes that the 
total intensity in every pixel is well modelled by a finite sum of 
separate emission mechanisms such that 



Itot (i', n) = ^ Ti (i/, h) 

i=l 
n 

= ^ 9i{u,n) Ti{voi,n) 



(2) 



(3) 



(8) 



If each analysis is limited to a single sky region covered by a single 
set of parameters 6i{y) then each of the individual terms on the 
right-handside of Equation l[8j are monopole terms. By extension, 
the total 5tot{v) must also be a constant monopole. Therefore, we 
can define a new template Toslv, n) = ToS = 1 at all frequencies 
and all sky positions. Equation Q can then be re-written 



(9) 
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where 6'(„+i) (u) = Stot (i^) and Tl^^+i) (") ^ ^off - The C-C analy- 
sis then proceeds just as in the case of a zero-offset but simply with 
an additional template term, and the value of Stot (i^) is returned 
along with the other n fit parameters, 9i{u). 

The introduction of an extra fitting parameter (the amplitude 
of the monopole-term) increases the uncertainties returned in the 
other parameters of the C-C analysis. This effect is the largest at 
high Galactic latitudes where there is little variation in the fore- 
ground emission. That is, one should expect the correlation of the 
data with the monopole-template to increase with increasing lati- 
tude even if a true constant-offset does not exist in the data maps. 

When these offsets are included in the C-C the reduced-x^ is 
significantly better than without the offset. At K-band it decreases 
by about 24% for the temperature analysis and by about 8% in the 
analysis of the polarised emission. These offsets are found to ex- 
hibit a distinct synchrotron-like frequency spectrum, whereas one 
would naively expect little correlation with frequency. This sug- 
gests that the Haslam et al. 408 MHz synchrotron template is dom- 
inating the fit offset terms. 



4.2 Uncertainties 

Writing the C-C error equation following the component formal- 
ism, the uncertainties on the fitted C-C coefficients for the i*'^ tem- 
plate are 



(10) 



In Equation JlOt . the covariance matrix Msn = Ms+Mn, where 
Mm and Ms are the covariances of the instrumental noise and the 
CMB, respectively. The CMB signal covariance only has a signif- 
icant contribution to the temperature analysis, where it dominates 
the instrumental noise. 

We compute the signal covariance matrix in the usual way, 
taking the CMB power spectrum Ci from the WMAP best-fit cold- 
dark-matter (ACDM) power-law spectrum model (Nolta et al. 
2009). In practice, we take into account spatial correlations within 
30° of a given pixel, effectively neglecting correlation values that 
are at most 0.5% that of the pixel auto-correlation. The noise co- 
variance in intensity is determined from the uncorrelated pixel 
noise as specified for each pixel in the WMAP data; this is an A'^ x 
diagonal matrix with the inverse of the observation number in each 
pixel on the diagonal (ctj — ao / \/Nobs)- 

We use the same strategy as in Eq. l llOt to compute the errors 
in the polarisation coefficients defined in Section 1431 (Eg . l23t . In 
the polarisation analysis the covariance matrix is dominated by the 
noise, which exceeds the contribution from the polarisation power 
spectrum Cf^ by two orders of magnitude.* Mn is, to a good 
approximation, a diagonal matrix as the pixel-to-pixel noise cor- 
relation is at most 1% of the diagonal term. We verified this by 
inspecting the 4 x N x N complete noise matrix available at the 
LAMBDA website, even though this is given at A^'side = 16 and 
our work is done using A^sido ~ 32 for the reasons introduced in 
Section |2] While it is true that we have no information on exact 
pixel-to-pixel correlations at A^sido ~ 32, the A^sidc ~ 16 noise 
matrix suggests there are no large-scale noise correlations for the 
polarisation analysis. 

The polarisation C-C results are, therefore, computationally 



On the diagonal elements, Ms = 0.45 /jK^. while Mn = 55 /iK^. 



easier to determine than the temperature results, because we only 
need to compute A'^-dimensional arrays as opposed to the N x N 
needed in the temperature case. Therefore, we compute the noise 
for the total polarisation in each pixel by combining noise in the 
Q and U maps at each frequency and properly taking into account 
Q-U correlations. 

In every pixel n we have a matrix 



A'(n) 



QQ{n) 
UQ{h) 



QU{n) 
UU(h) 



where QQ[n) and UU{h) are noises and QU(n) is the n pixel 
Q-U noise correlation, all in terms of observation number. Uncer- 
tainties on Q and U are given by the inverse of N{n) 



Propagating these noise terms to compute the uncertainty on P = 
yCP+lj^ we obtain: 



(jp{n) 



(11) 



The (Tp(n)^ are the elements of the (diagonal) noise matrix Mn 
for the total polarisation P. The uncertainty on the polarised corre- 
lation coefficients is: 



5pi{v) 



Y,[{TT{h) ■ m^'\v, h) ■ T{h))-^]u. (12) 



When choosing the area over which to perform the C-C analy- 
sis, one should consider the trade-off between the uncertainties and 
the accuracy of the modelling. While it is more likely that in a small 
area a single fitting coefficient is sufficient to represent a given fore- 
ground, the uncertainties calculated with Eq. l llOt and M2\ decrease 
with the square root of the number of pixels considered. In the fol- 
lowing, our choice to consider relatively large areas is a trade-off 
between these two competing aspects. 



4.3 Polarisation analysis and fractional polarisation 

When dealing with the cross-correlation of the polarised emis- 
sion some preliminary comments are necessary. Since Equation Q 
holds for any measured intensity and polarisation value, if tem- 
plates from all Stokes parameters were available at some set of 
reference frequencies, i.e., Qi{iyoi,h) and Ui{yoi,n), then Equa- 
tion l[4j could be solved separately to determine the coefficients 
Oiiy) in polarisation. 

That is. Equation l|4]l, written for Q and U in each pixel n 
becomes 



Qtot(t^, = Qi{v,n) 
1=1 

n 

Utot{v,h) = ^ Ui{iy,n) 



(13) 
(14) 



Inserting the definitions of the Stokes parameters in terms of the 
fractional linear polarisation / and the polarisation angle Q — 
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fl cos(27) and U — fl sin(27). Equations J13I ) and J 14b become 

n 

Qtot{v,h) = ^ qi{v,n) Tt{v,n) (15) 



Utot{i',n) = ^ Ui{u,n) Ti{v,n) 



(16) 



where the correlation coefficients are 

gi(!/, n) = 7i)6li(!/) cos[27i(i^, 77,)] and (17) 
Ui{v,n) ^ fi{i/,h)9i{i/) sin[27i(!/, n)]. (18) 

To apply C-C as we did in the temperature analysis, we assume 
again that the coefficients qiiy^h) and Ui{u,n) are constant over 
the analysed region of the sky. 

We could still apply this procedure if we had complete in- 
formation about the polarisation angle for every foreground. For 
example, one could use simulated Galactic magnetic-field maps to 
estimate the position angle of a thermal-dust polarisation-template 
and the K-band polarisation data as a synchrotron angle template 
as has been done by other authors (e.g., Dunkley et al. 2009; Page 
et al. 2007). 

In this work we wish to determine if useful knowledge can 
be extracted from the WMAP polarisation data without making any 
assumptions about the dust or synchrotron polarisation angles. It is 
possible to do so by considering the total polarisation 

PL{y,n) = f{u,n)T\u,n) = Q\u,n) + U\p,n) (19) 



= Q^g,(//,n)r.(n) + [Y,u,{u,n)T,{n 



(20) 

Equation |20] can be re-written by replacing qi and m using 
Equations | |17I > and so that 

n 

n n 

+ E E hfje^e,T,{n)T,{h) cos [2 (7,(n) - 7,(n))] . (21) 

(Note that we have dropped the frequency-dependence in Equa- 
tion 1211 for simplicity only.) If we further assume that the polarisa- 
tion angles are equal for all emission components (7i(n) = 7j (n)) 
then 



where 



Ptot{y,n) ^^p,{u)T,{f, 



Pi(u) = fi{i')ei{v). 



(22) 



(23) 



Equation l l22b is analogous to Equation but with the 
polarised-intensity correlation coefficients rather than 6i[v). 
Therefore, Pi(i^) can be computed by analogy with the C-C used to 
calculate the 9i{v) given by Equation (|5}. The fractional polarisa- 
tion coefficient fi(p) can be extracted for each component at each 
frequency from these two correlation coefficients: 



with an uncertainty of: 



5pi 



+ 



(24) 



(25) 



The assumption of an equal polarisation angles for all emis- 
sion components is reasonable because the direction of all dif- 
fuse emission types is determined by the direction of the mag- 
netic field. This assumption leads to a total polarisation (summed 
over all emission components) with angles equal to the individ- 
ual components, so that it is easily falsifiable by any measurement 
in which the total position angle changes with frequency. While 
such angle rotations with frequency are observed in the WMAP data 
they are typically smaller than ~ 20° at a resolution of 4° . More- 
over, the nearly-constant angle assumption is supported by Kogut 
et al. (2007) who found that the fits based on K-band (synchrotron- 
dominated) polarisation angles or the starlight-polarisation angles 
gave similar results, with a difference in the goodness-of-fit of 
about 1%. The advantage of our approach is that no models are 
required to estimate the polarisation angles. 

The fractional polarisation fi(i^), as given in Equation i24l . is 
the ratio of correlation coefficients. One should keep in mind that 
negative values of the correlation coefficients pii^v), and therefore 
fi{i^), are in fact physically meaningful. Anti-correlations between 
the polarisation data and a total intensity template will occur in 
regions where the polarised emission drops with increasing total 
intensity. For example, this behaviour can be expected in regions 
where the polarisation drops due to increased line-of-sight averag- 
ing towards regions of increasing column density. This point pro- 
vides another reason to consider small areas of the sky where the 
polarisation fraction is more likely to be constant. 



4.4 Polarisation noise bias 

Given the Q and U values in each pixel, we calculate the polarised 
intensity using 



P(n) = ((Q)^(n) + {Uf(n)y^\ 



(26) 



and its associated uncertainty, where the direction h denotes each 
A'sidc ~ 32 (~ 2° X 2°) piece of sky. The correlation analysis is 
then carried out on the P(n) in each section of the sky defined in 
Section|5] 

The polarisation is a complex quantity which can be described 
by either its real and imaginary parts (Q and U), or its amplitude 
and phase (P and 7). If one assumes gaussian statistics for Q and 
U, then the non-linear transformation into P and 7 results in non- 
gaussian statistics for these quantities. The key point for our analy- 
sis is that P has an asymmetric probability distribution (e.g., Vail- 
lancourt 2006; Simmons & Stewart 1985) resulting in a positive 
bias. We have corrected for this bias using the formula: 



ftruc — 
-Ptruc ^ 



for P/ap < V2 (27) 
P2 - cr2 for P/ap > V2 (28) 



where ap is the polarisation uncertainty defined in Section \4m 
While the last expression is not exactly precise one can limit the ef- 
fects of the bias by only using data above some polarisation signal- 
to-noise level. ^ Such a bias can be estimated using monte carlo 
simulations as described in Section lSTI 



No bias-correcting algorithm is exact. For more information see Sim- 
mons & Stewart (1985), Vaillancourt (2006), and references therein. 
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5 RESULTS 

We used C-C to compute, at each WMAP band, the emission coef- 
ficients of the three Galactic foreground components plus the offset 
monopole template. The analysis presented in this work was done 
in three different spatial configurations: on the whole sky outside 
of the KQ85 Galactic mask (hereafter the "all-sky" analysis), in 
different isolatitude slices (hereafter the "latitude" analysis) with 

6 = -907 - 50°, -507 - 20°, -20°/20°, 20°/50°, 50°/90°, 
and in 48 HEALPix A'^sidc ~ 2 regions (hereafter the "pixel" analy- 
sis). Selecting such partitions allows us to investigate how the fore- 
ground properties depend on direction and latitude (see also Kogut 
et al. 2007; Miville-Deschenes et al. 2008), while increasing the 
likelihood that the area we are analysing has a uniform fractional 
polarisation (see Section l43t . The best pixel area size among those 
tried is A'sidc = 2; it is small enough that the assumption of con- 
stant polarisation fraction is more likely to be satisfied, and large 
enough to contain spatial variations needed for the C-C method. 
Moreover, these areas match with previous works (e.g. Dunkley et 
al. 2009) allowing us to compare our results. Fig.[T]shows the sky- 
locations of pixels used for the pixel-analysis and the latitude slices 
used for the latitude-analysis (see Table [3] for the mean latitudes 
and longitudes of each pixel). 

From the C-C coefficients we derive results on the relative 
level of each foreground in a given region, its frequency depen- 
dence in temperature and polarisation, and its polarisation fraction. 
We note here that it is not always possible to extract all this infor- 
mation at all frequencies and for all regions due to limited signal- 
to-noise ratios. Therefore, we concentrate on results where there is 
2(T statistical significance. 

5.1 Tests on a simulated sky 

In order to test the analysis method described in the preceding sec- 
tions, we ran the CCA code on 500 noise realisations of a simu- 
lated sky with a given foreground parametrisation. We created sky 
maps in both temperature and polarisation, based on the three tem- 
plate temperature maps from Section [3] and scaled them using the 
average coefficients from Davies et al. (2006) at each frequency. 
We used a constant polarisation fraction of 20% for synchrotron, 
3.5% for dust^, and 0% for free-free and a polarisation angle for 
both synchrotron and dust based on the WMAP5 K-band map. We 
also made a set of simulations with no dust polarisation to compare 
cases without and without dust polarisation. The CMB model is a 
Gaussian realisation based on the WMAP5 concordance model. We 
used the LAMBDA A'sidc = 512 /, Q, and U inverse covariance 
matrices to produce a A'sidc = 512 noise realisation, which was 
added to the CMB and foreground maps, and then degraded to the 
resolution used in our analysis, i.e.. Aside = 32. 

We compared the results for each of the three components for 
the 500 realisations in each pixel. We focused on the results from 
K-band, where the data had the most significance. The recovered 
coefficients were, in general, close to the input values with no large 
biases. In temperature, the results were all within 3a of the input 
value, except for one region for synchrotron (region 10 at 3.3(t) and 
one region for dust (region 12 at 3.1a). The estimated uncertainties 
corresponded well to the scatter in the coefficients from the 500 
realisations. 

^ The 3.5% polarisation represents the polarisation for the total dust- 
correlated component at K-band, which will be dominated by anomalous 
microwave emission, and virtually no contribution from thermal dust. 



In polarisation, the results for synchrotron and free-free were 
in good agreement with the input values. The mean recovered frac- 
tional polarisation for synchrotron was 20.0 ± 0.3 % for the case 
with no dust polarisation and 19.5 ± 0.3 % with dust polarisation. 
For free-free, we obtained mean values of —0.02 ± 0.18% and 
—0.04 ± 0.19% for each case, respectively. This gives us con- 
fidence that the method works as expected and yields reasonable 
results. 

However, the values recovered for dust did indicate bias i.e. an 
additional systematic error. The recovered fractional polarisation 
was 0.003 ± 0.15 % and 0.38 ± 0.15 % for the simulations with 
no dust polarisation, and 3.5 % dust polarisation, respectively. This 
shows that simply averaging all the polarisation values within a 
pixel will not yield a reliable estimate if the true value is as low 
as a few percent. The situation improves when only using regions 
with significant (> 2a) detections. We found 9 regions that were 
significant at the 2a level, which resulted in a weighted average 
of 2.0 ± 0.2 %, i.e., an underestimate by 1.5 % in absolute value. 
For the case with no dust polarisation, no pixels were detected at 
greater than the 2a level. Similarly, for free-free, we obtained an 
average value of 1.4 ± 0.4 % when using the 6 regions that were 
significant, i.e., an overestimate of 1.4 % in absolute value. 

From these tests, we conclude there is a bias in the average 
fraction polarisation of ~ 1.5 %, when the true fractional polarisa- 
tion is a few percent or less. For a single region, we found cases 
where the bias is up to ~ 5 %. The reason for this bias appears 
to lie in cross-talk between the various components. For example, 
for regions we recover a significant dust polarisation, with no dust 
polarisation in the simulation, one can see a bias in either the syn- 
chrotron or free-free polarisation; this can explain why the average 
synchrotron polarisation is slightly low, at the 1.7a level, in the 
case where there is dust polarisation. Nevertheless, we can still ob- 
tain useful constraints on the fractional polarisation with the current 
data, keeping in mind the amplitude of possible biases. 

Also, we can use the results of these simulations to estimate 
the significance of results for individual pixels. These will be dis- 
cussed separately, in the following sections on the synchrotron, 
free-free and dust results. For the cases where there is significant 
bias observed in the simulations, the fractional polarisation results 
will be presented in terms of a statistical (stat) and systematic 
(sys) error. For example, for the significant dust regions above, the 
weighted average is 2.0 ± 0.2 (stat) ± 1.5 (sys) %. 

5.2 All-sky versus pixel analysis 

We show in Fig. [2] the contribution of each foreground to the tem- 
perature emission as a function of frequency. These estimates, re- 
ported in Table [T] were obtained by multiplying the C-C coeffi- 
cient for a given foreground at a given frequency with the root- 
mean-square (r.m.s) of each foreground template outside the KQ85 
mask. The figure clearly shows that the dust-correlated emission is 
the dominant component in temperature at all frequencies (see also 
Davies et al. 2006). The rise from 61-to-94 GHz indicates the high 
frequency contribution of thermal dust, while at lower WMAP fre- 
quencies, there is a significant contribution from anomalous dust. 
It is this low frequency dust-correlated emission that has been ten- 
tatively identified with emission from small rapidly-spinning dust 
grains (Draine & Lazarian 1998a,b). However, the dust contribu- 
tion to polarisation is typically subdominant at low frequencies and 
will be discussed further in Section [5!4l 

Given the C-C coefficients and foreground templates, and the 
assumption that the templates trace all the foreground emission, we 
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Figure 3. Residual maps (on mollweide projection) for the WMAP K-band temperature (top) and polarised intensity {bottom) analyses. The results using the 
all-sky (KQ85 cut) region are reported (left) together with the results using the 48 HEALPix A^side = 2 regions (right). The units of the color-scale are /^K. 




Figure 1. All-sky map with the 48 HEALPix N^idc = 2 regions in dif- 
ferent colours with the KQ85 mask applied (dark blue). Horizontal black 
lines represent the isolatitude slices (latitude analysis) with b = —90°/ — 
50°, -50°/ -20°, -20°/ 20°, 20°/ 50°, 50°/ 90°. 

can generate a model of the total foreground emission in the WMAP 
data. Comparing this model to the actual WMAP data yields a set of 
residuals which we use to estimate the extent to which a single fit- 
coefficient (in the case of the all-sky analysis) represents the fore- 
grounds across the whole sky. Results for the K-band, where the 
foregrounds are the brightest, are reported in Fig. [3j/e/i') for both 
temperature and polarisation. In the all-sky temperature residuals, 
the largest residuals occur near the Galactic plane, particularly near 
the Galactic centre. This feature has been interpreted as the Galac- 
tic haze, possibly associated with annihilating particles (Finkbeiner 
et al. 2004; Cumberbatch et al. 2009). However, this feature is not 



100 




20 30 40 50 60 70 80 90 100 

Frequeny (GHz) 



Figure 2. R.M.S. emission spectrum of Galactic foreground components 
from the all-sky (KQ85 cut) intensity and polarisation analysis, in /iK. Only 
the polarised synchrotron emission is plotted as the dust-correlated and free- 
free-con'elated polarised emission are consistent with zero for the all-sky 
region. 

evident in the pixel-analysis where we compute separate C-C coef- 
ficients in 48 different regions (P'lg^right)). This is a consequence 
of a varying synchrotron spectral index, which is somewhat flatter 
than average near the Galactic centre, due to a diversified popula- 
tion of electrons in the Galaxy. 

Bright residuals in the all-sky analysis are also concentrated 
on the well-known Galactic spurs such as the North-Polar Spur 
(NFS) above and below the Galactic centre region. These spurs are 
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K-band 


Ka-band Q-band V-band 


W-band 


Temperature r.m.s. emission [fiK] 


S 
D 

F-F 


54.8 ±»- ' 
75.0 ±0-6 
32.0 ±0-5 


15.8 ±"- ^ 7.2±»-« 0.8 ±1 
27.6 ±o « 16.1 ±0-7 10.1 ±0-9 
14.6 ±0-5 8.5 ±0-5 2.0±0-e 


-0.8 ±1=' 
13.6 ±1-3 
-0.2 ±1-1 


Polarisation r.m.s emission [^iK] 


S 
D 

F-F 


11.9 ±"-2 
-0.2 ±"1 
-0.2 ±01 


3 9 ±0.2 J 7 ±0.3 1 ±'*-'' 
-0.1 ±0-2 o.O±0-» 0.1 ±0-3 
-0.3 ±0-2 -0.2 ±0-2 0.0 ±0-4 


0.1 ±"-4 

0.0 ±0 * 
0.0 ±0.6 



Table 1. Emission r.m.s. values for the all-sky analysis in both temperature 
and polarisation for synchi'otron (S), dust (D), and free-free (F-F) in the five 
WMAP bands. These values are plotted in Fig.|2] 





K-band 


Ka-band 


All-Sky Temperature 


0.94 


0.68 


All-Sky Polarisation 


5.37 


0.73 


Latitude Temperature 


0.97 


0.74 


Latitude Polarisation 


3.72 


0.56 


Pixels Temperature 


0.65 


0.48 


Pixels Polarisation 


1.65 


0.30 



Table 2. Reduced-x2 for the three analysis types implemented in the paper. 
The values are the average over the different regions (when more than one 
is used). 

known to have a slightly steeper synchrotron spectral index than 
average and, therefore, are not well suited for the all-sky polari- 
sation analysis. Furthermore, these large coherent features are ex- 
pected to have significantly different fractional polarisations. The 
residual map for the 48 pixel-analysis (Fig.O shows less extended 
structure, with more localised residuals near the mask edge. This 
indicates that a local analysis is more suited to handle the complex 
foreground variations. In summary, Fig. [3] shows that it is subop- 
timal to model foregrounds without considering that they are spa- 
tially varying in the sky. 

To have a quantitative understanding of the goodness-of-fit we 
compute the reduced-x^ for each type of analysis, both in temper- 
ature and polarisation. These are reported for the WMAP K- and 
Ka-bands in Table |2] Again we see that neither the all-sky nor the 
latitude analyses are as accurate as the pixel-analysis, especially in 
the K-band polarisation. A limitation of the pixel-analysis is that, 
in any given pixel, there may not be sufficient data points to yield 
a signal-to-noise ratio larger than 1. This is a particular problem 
for high frequencies (mostly V- and W-bands) in polarisation. In 
most cases we will be limited to studying the WMAP foregrounds 
at K-band only. The 48 emission coefficients in K-band, for each 
foreground component in intensity and polarisation, are reported in 
Tabled 

5.3 Synchrotron 

Synchrotron is the dominant contribution to the polarised K- and 
Ka-band emission, especially in the northern hemisphere, as re- 
ported in Table [3] because of the strong emission associated with 
the north polar spur. On the other hand, synchrotron is comparable 
to anomalous dust emission in temperature, and again its emission 
is more intense in the northen part of the sky. 

From the all-sky analysis, using K- and Ka- band and the cor- 
rected Haslam et al. map described before, we find /?s = —3.32 ± 
0.12 in intensity and /3s = —3.01 ± 0.03 in polarisation, where 
the errors come from the uncertainties on the C-C coefficients. The 
smaller polarisation uncertainties occur here because synchrotron 



is the dominant polarised foreground at low frequencies, as op- 
posed to the intensity analysis where the dust-correlated and free- 
free-correlated emission are comparable to the synchrotron. These 
synchrotron spectral indexes are in good agreement with previous 
results, i.e. /3s ~ —3 (e.g. Miville-Deschenes et al. 2008). 

From the pixel analysis we find that the synchrotron emission 
is polarised from « 5% in some regions at low Galactic latitude 
up to ~ 40% at high latitudes, as reported in Fig.|4](/cp). In Fig.|5] 
the same quantity is plotted as a function of the mean Galactic lat- 
itude of each region. The decline of the polarised emission at low 
latitudes is well known in the literature (e.g., WoUeben et al. 2006) 
and is interpreted as a depolarisation effect near the Galactic plane 
(|6| < 20°) due to integration along the line-of-sight of emission 
with different polarisation angles. The average polarisation fraction 
in the pixel regions with ( 1 6| < 20°) is 8.6±1.7 (stat)±0.5 (sys) %. 
Uncertainties at high latitudes are usually larger than at low latitude 
because the templates have less variations in their structure at high 
latitudes. This renders the C-C analysis increasingly degenerate, as 
explained in Section l42l Nevertheless, the synchrotron fractional 
polarisation at high Galactic latitude is in the range ~ 10 — 40% 
with an average value of 19.3 ± 0.8 (stat) ± 0.5 (sys) % and stan- 
dard deviation of ±11.0% (|61 > 20°). As discussed in SectionlSH 
the synchrotron values are the most robust of the three components, 
and we believe the errors from the C-C method are representative, 
even after averaging. 

Our estimate of the synchrotron fractional polarisation is con- 
sistent with other work at the la level. Kogut et al. (2007) com- 
puted the fractional synchrotron polarisation by dividing the po- 
larised emission (estimated using a pixel-by-pixel frequency fit), by 
a synchrotron intensity map (based on a MEM analysis; Hinshaw et 
al. 2007). They find a typical fractional synchrotron polarisation of 
5-25 %. However, since this intensity map was computed ignoring 
a possible contribution from anomalous dust, the synchrotron in- 
tensity may be overestimated, resulting in an underestimate of the 
fractional polarisation. On the other hand, Miville-Deschenes et al. 
(2008) included a spinning-dust contribution to the total K-band 
intensity, and obtained a higher synchrotron fractional polarisation, 
up to 40%. The fact that our results are compatible with these im- 
plies that in most of the pixel regions the assumptions described in 
Section l431 are reasonable. 

We have computed the synchrotron spectral index in every 
pixel with a S/N ratio > 2 again using only the K- and Ka-bands; 
see Table |4l The problem in this case is that the pixels with S/N 
ratio > 2 are few, and we conclude that the C-C technique is not 
appropriate for a detailed spectral-index analysis. Nonetheless, the 
spectral indices in polarisation are consistent with other analyses 
(e.g., Dunkley et al. 2009) who find /3s — —3; we find an average 
value is /3s = -3.24 ± 0.20. 

In temperature, the results are not compatible with polarisa- 
tion, with an average spectral index of /3s — —2.02 ± 0.20. This is 
clearly too flat and must be an artifact of the analysis. We checked 
that these results, especially in the 5 low /3s regions, are not due 
to a high correlation between synchrotron and dust. Without suf- 
ficient differences in spatial morphology between synchrotron and 
dust, confusion can occur between the two components in the C-C. 
As we expected, the dust and synchrotron templates are correlated 
with an all-sky (outside the KQ85 mask) correlation coefficient of 
r — 0.48. In the regions where high synchrotron spectral indices 
are observed (4, 5, 16, 27 and 32) the correlation between dust and 
synchrotron is higher (r = 0.71). This suggests again that, in tem- 
perature, the C-C in some pixel regions is not able to correctly 
separate the synchrotron, free-free and dust components. In fact, 
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Pixel 


Pixels 


(Lat) 


( Long ) 


Synch K 


DustK 


Dust W 


Synch % 


Dust % 


F-F % 


Xt 


xi 


1 


207 


65.7 


46.2 


38.9 ±3.« 


19.6 ±4-^ 


11 ±5.3 


40.6 ±4 7 


< 11.0 




0.45 


0.41 


2 


214 


66.7 


134.8 














0.31 


0.81 


3 


214 


66.8 


224.5 




9.5 ±3.2 










0.49 


0.73 


4 


202 


67 


318 


30.6 ±3 * 


9.5 ±3-4 




23.8 ±3-^ 




< 11.6 


0.84 


0.41 


5 


189 


42.3 


21.3 


49 5 ^'^-'-^ 


56.8 ±4.3 




26.1 ±^-4 


7.2 ±^-3 




0.45 


3.01 


g 


195 


43 


69.7 




19.9 ±3-2 


6.2 ±3-5 




10.6 ±4 9 




0.49 


2.11 


7 


209 


42 


110.3 


9.9 ±2.6 


32 7 ±3.4 


7.4 ±3.7 


17.8 ±7-5 


< 5.0 




0.55 


0.75 


8 


195 


41.9 


159.8 


10.2 ±4-5 


30.8 ±4-3 


8.4 ±4-8 


42.4 ±21.6 


< 7.6 




0.49 


1.42 


9 


211 


41.7 


200.1 


8.7 ±3-7 






16.7 ±14-4 






0.56 


0.22 


10 


232 


41.6 


250.8 


12.9 ±3.4 


12.8 ±2-8 




18.3 ±7 '' 


< 12.8 


< 12.4 


0.60 


0.31 


11 


205 


42.1 


289.8 


16.8 ±2.9 


23.8 ±3-4 




21.3 ±S-2 




< 11.6 


0.47 


1.05 


12 


214 


42.4 


339.6 


25.7 i"*-! 


43.5 ±4-2 




22.5 ±5''' 


< 5.8 




0.37 


2.16 


13 


44 


30.9 


4.1 


21.1 ±4-5 


69.6 ±4 9 


14.1 ±5-2 


44.3 ±12.8 


< 3.6 


< 6.4 


1.32 


2.32 


14 


200 


21.1 


44.9 


104.8 


54.7 ±"'•3 




23.7 ±1.4 






0.98 


2.02 


15 


158 


21.4 


88.6 




110.5 ±3 4 


17.2 ±3-4 




< 1.6 


< 5.0 


0.63 


1.87 


16 


167 


22.7 


135 


81.9 ±5 3 


76.4 ±4-2 




19.8 ±19 




7.6 ±4-8 


0.61 


0.28 


17 


211 


20.1 


180.5 
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48 
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313.6 
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Table 3. Results of the pixel analysis. Here we report the pixel region number, the number of data points, the average latitude and the average longitude of 
each region, the synchrotron-correlated temperature emission at K-band in /xK, the dust-coiTclated temperature emission at K-band and at W-band in /.jK, and 
fractional polarisation percentage of synchrotron, anomalous dust and free-free at K-band. Bold faces are used to emphasise a component with a temperature 
emission higher than 60 at K-band and 15 ^K at W-band, and are use to indicate regions where a particular component is found to be more than 2a 
significantly polarised. We indicate 2a upper limits when in that region the fractional polarisation is compatible with zero. No data for fractional polarisations 
indicates upper limits reaching 100%. In column 5, 6 and 7 we display detections that are significant at the la level. Note that the systematic uncertainties are 
not included here and can be significant when the polarisation fractions are low, i.e. less than a few per cent (see Section lsTTI . The values from temperature 
and polarisation analysis at K-band are also reported. 



most of the areas where the spectral index is biased by this effect 
are the ones that have problems for synchrotron in our simulations, 
e.g. regions 16, 27, 32. Moreover, when fitting over larger regions, 
the spectral index moves to more typical values of /?s ~ —3. The 
same difficulty with C-C analysis was found in previous work, e.g. 
Davies et al. (2006), where the introduction of the 0.4 GHz point 
of Haslam et al. was needed to obtain meaningful spectral index 



values. However, we note that this problem is most significant at 
the higher frequency channels above 33 GHz while K-band the co- 
efficients are not strongly affected. Indeed, the spectral index com- 
puted from the weighted mean K-band coefficient alone, which cor- 
responds to the mean spectral index between 408 MHz and 23 GHz, 
is/3 = -2.97 ±0.01. 
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Table 4. Synchrotron spectral index values in different pixel regions from 
the intensity and polarisation analyses at K- and Ka-bands. Only the regions 
with a significant result are reported (cr^g < 1.0). 

5.4 Dust 

It is clear from the all-sky analysis of Fig. [2] that the total un- 
polarised signal across all 5 WMAP bands is dominated by dust- 
coirelated emission. The rise in dust emission from 61 to 94 GHz is 
consistent with the expectation of thermally emitting dust. The drop 
with frequency from 22 to 61 GHz has been termed "anomalous" 
dust-correlated emission and has been attributed to such sources 
as spinning-dust emission (e.g., Lazarian & Draine 1998a,b) and 
dust-correlated synchrotron emission (Bennett et al. 2003). The 
results in Fig. |2] and K-Ka spectral index for anomalous dust of 
/3(j ~ —2.5, is consistent with the results of Davies et al. (2006). 

The results for emission presented in Table [3] show that dust 
emission also dominates the unpolarised K-band emission for most 
of the sky in the pixel-analysis. We find that the dust-correlated 
emission is dominant at low frequencies in many pixels even where 
the thermal dust emission is not present, (e.g., pixels 16, 22, 23, 27, 
and 40). 

Fig.|4]shows the polarisation fraction of each foreground at K- 
band, computed using Equation i24l . The same values are reported 
in Table [3] As expected, synchrotron emission is highly polarised, 
while dust and free-free almost unpolarised. We compared the val- 
ues found in various pixels with the distributions obtained from 
simulations. As a general comment on the simulations we can say 
that a) the range of the recovered dust polarisation fraction, assum- 
ing the null hypothesis, is typically 5-10% and varies between a 
few percent and 25% (for regions very distant from the Galactic 
plane, e.g., 2, 3, 47, 48); b) the distributions from simulations with 
and without dust polarisation are, in general, very similar, although 
higher mean polarisation values are not always associated with the 
case corresponding to the polarised dust (this may indicate cross- 
talk with other components); c) the variance in the distributions is 
typically similar to the error bars computed by the CCA on real 
data. 

When we compare the detected values of the polarisation frac- 
tion in the K-band with the distribution we find that in some pix- 
els they exceed the tail of the simulation distribution. The pixels 
for which the simulation-recovered values are always below the 
(real) detected dust polarisation are 5, 17, 21, 23, 25, 27, 34, 35 
(see examples in Fig. |7]l^. In some regions (5, 6, 17, 21, 23, 24, 

^ For comparison, the number of pixels in which the fractional polarisation 
is always below the distribution is two. 





Figure 4. Synchrotron-correlated (top), dust-correlated {middle), and free- 
free-correlated (bottom) fractional polarisations in the 48 regions of the K- 
band pixel-analysis. The regions compatible with a zero polarisation frac- 
tion within 1(T are plotted in blue. Regions with 2ap > 30%, or with neg- 
ative fractional polarisation, are plotted in grey. 

25, 27, 32, 34, 35, 45) the detections are at more than the 2a 
level. Some of these regions also show a (small) increase in syn- 
chrotron polarisation with respect to simulation results (e.g., 5, 17, 
32), with an accompanying anti-correlation with the dust values. 
Indeed, in some cases, the recovered dust polarisation fraction be- 
comes slightly negative, such as in region 5 (see Fig. IT}, which is 
evidence of cross-talk between components. Six regions (5, 17, 23, 
34, 35, and 45) out of the 48 exceed the simulated distributions at 
the 99th percentile and are individually significant at > 2a. They 
showed fractional polarisations in the range 2.6-7.2 % except for 
the anomalous pixel 45, which appears to have a high polarisation 
fraction of 32 ± 12%, detected at a significance level of 2.7ct. This 
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Figure 5. Synchrotron- and dust-correlated fractional polarisation of 
WMAP K-band as a function of mean Galactic latitude in the 48 regions 
of the pixel-analysis. We plot only regions where /p > crjp . 



c 
o 



E 
tu 

o 



JO 

o 

Q. 



Pixel 5 
Pixel 21 
Pixel 23 
Pixel 25 
Pixel 27 
Pixel 34 




40 50 
Frequency (GHz) 

Figure 6. Dust-coirelated polarised emission r.m.s. in as a function of 
frequency in the pixel-analysis. We report only the six pixel regions where 
there is a 2cr significant emission in the K-band. 



region does not contain any strong foreground features in total- 
intensity or in polarisation. Additionally, the simulation recovers 
the correct polarisation but with a relatively large uncertainty. The 
detected value in the real data is well outside the distribution, as 
shown in Fig. [8] This region should be investigated further with 
higher sensitivity data (e.g., WMAP 9-year, Planck). 

Although we cannot rely heavily on individual pixels, we can 
still continue, knowing that there can be biases of up to ~ 5 % in 
any given region, and ~ 1.5 % when averaging significant (> 2(t) 
pixels, as discussed in Section |5T| Furthermore, we can choose 
pixels that in the simulations were able to recover the correct dust 
amplitude within la of the true value (regions 6, 12, 25, 33, 34, 35, 
39, and 45). Although these may not be completely representative 
of the real data, they should be more reliable than simply averaging 
all or just significant pixels. For these, we find a weighted average 
of 3.2 ± 0.9 (stat) ± 1.5 (sys) %. We do not claim this as a real de- 
tection, since this is only a ~ 2(j detection when considering the 
overall bias seen in the simulations described in Section lSTI How- 
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Figure 8. Histogram of the K-band dust polarisation fraction for region 
45, from 500 realisations, for the case with no dust polarisation (red line) 
and with 3.5 % dust polarisation (blue line). The input values are shown as 
vertical dotted lines. The green line and error-bar represent the best-fitting 
value and its la uncertainty from the real data. 



ever, it is clear that the average K-band dust polarisation fraction is 
less than about 5 % at 95 % confidence. 

A comparison of the dust and synchrotron values at K-band 
(Fig. [5j as a function of Galactic latitude) shows that the K-band 
synchrotron polarisation increases from fs ~ 5% to ~ 40%, while 
the dust remains fairly constant within the error-bars at fd ~ 2% — 
10% in the pixels where significant dust polarisation emission is 
found. Over most of the sky dust is found to be unpolarised, with a 
2(T upper limit on fractional polarisation of about 5% in seventeen 
pixel regions and about 10% in four other regions. 

Given the limitations of our results, particularly at frequen- 
cies above 33 GHz, we cannot make a detailed analysis of the fre- 
quency dependence of the dust polarisation. Fig.|6]reports the r.m.s. 
of the dust polarisation emission. The fractional polarisation for 
these same data are plotted in Fig. |9] (pixel 21 has been removed 
from this plot as the temperature-analysis yields negative C-C co- 
efficients). Here we report just K-, Ka,- and Q-bands, as the higher 
frequency data have signal-to-noise ratios less than la. While the 
fractional polarisation data do exhibit a slight frequency depen- 
dence the uncertainties are quite large. Due to these large uncer- 
tainties we provide only a broad discussion, rather than a detailed 
analysis, of these results in the remainder of this section. 

In the all-sky analysis (Fig. |2j the total emission from the 
anomalous-dust component appears to dominate the thermal dust 
component for frequencies below 61 GHz (V-band). We find simi- 
lar results in the temperature data of the five pixels of Fig. |9] (not 
shown). Lazarian & Draine (2000) predicted that spinning-dust po- 
larisation should be no greater than 2% at 20 GHz and become al- 
most completely unpolarised above 40 GHz. If the anomalous com- 
ponent is due to spinning-dust then the dust-correlated polarisation 
should be consistent with zero-polarisation in the frequency range 
40-61 GHz (Q-V bands), with the possibility of increasing polar- 
isation towards lower frequencies. Within the error-bars, this is at 
least consistent with the behaviour observed in Fig.|9] Alternatively, 
Draine & Lazarian (1999) predict polarisations as high as 40% in 
the 1-200 GHz range for emission from vibrating magnetic-grains 
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Figure 7. Histograms of the K-band dust polarisation fraction for tliree regions, from 500 realisations, for the case with no dust polarisation (red line) and with 
3.5 % dust polarisation (blue line). The input values are shown as vertical dotted lines. The green line and error-bar represent the best-fitting value and its Icr 
uncertainty from the real data. 
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Figure 9. Dust-coirelated polarisation fraction as a function of frequency 
in the pixel-analysis. We report only the five pixel regions where there is 
a > Icr result in the K-band that also have a non-zero value at Ka and Q 
bands. 



percent, apart from the regions around the Galactic poles that are 
essentially unconstrained. We find that in region 16 the detected 
free-free polarisation is higher than all simulation results but still 
compatible with zero at the 2cr level when C-C errors are consid- 
ered. In 16 of the 48 pixel regions we were able to put 2(7 upper 
limits on fractional polarisations ranging from 10% (in seven re- 
gions) to about 5% (in nine regions). These upper limits are more 
stringent where the uncertainties are small, i.e., where the free-free 
emission is higher, such as near the Galactic plane at longitudes be- 
tween 200° and 300°. Some pixels (e.g., 14, 16, 17, 24, 31, 33, 40, 
and 43) do show significant cross-talk in the simulations. By com- 
paring the simulations with the real data, we can isolate regions 
that we believe are not strongly affected by cross-talk and where 
the C-C errors are comparable or larger to the scatter in the simu- 
lations (therefore no significant bias). These are 4, 10, 11, 13,15, 
19, 22, 23, 26, 27, 32, 34, and 41. Examples are shown in Fig.fTOl 
The weighted average of these more reliable pixels is 0.0 ± 1.7 %, 
corresponding to an upper limit of 3.4 % (95 % confidence). The 
standard deviation in all 48 pixels is ±5.5%. We note that, unlike 
for dust, no regions exceed the simulation distributions at the 99th 
percentile and are individually significant (at > 2cr). 



(although the actual polarisation is not expected to be so high given 
some level of beam depolarisation and depolarisation due to line- 
of-sight integration). Depending on the exact number of magnetic- 
domains and the grains' aspect ratios, the magnetic-grains' polari- 
sation may drop from its maximal value to zero almost anywhere 
within the same frequency range; the observed polarisation as a 
function of frequency could have almost any behaviour (i.e., rising, 
falling, flat, or all three). Given the large uncertainties, we can only 
note here that the data plotted in Fig.|9]is not inconsistent with this 
behaviour. 

5.5 Free Free 

A simple fit to a power-law spectrum for the free-free emission 
in the all-sky analysis yields a free-free spectral index of /?/ = 
—2.15 ± 0.22. In the latitude analysis and in the all-sky analysis 
we find free-free emission is unpolarised (< 1%) over the sky. This 
result is confirmed in the pixel-analysis, where over most of the sky 
free-free is unpolarised at the Icr level, as reported in Table[3] 

The upper limits in the distributions of the polarisation frac- 
tion from simulations in K-band typically range between 5 and 15 



6 CONCLUSIONS 

We computed the cross-correlation of the intensity and polarisa- 
tion WMAP data in different sky regions using templates for syn- 
chrotron, dust, and free-free emission, plus a monopole-offset tem- 
plate. This technique is reliable because it takes into account the 
correlations between foregrounds, and because it requires no as- 
sumptions about the foregrounds' frequency behaviours. However, 
this study relies on the assumption of equal polarisation angles for 
the three foregrounds and on the hypothesis of constant fractional 
polarisation fraction in each of the the analysed sky-regions. We de- 
rive the frequency dependence and polarisation fraction for all three 
components in 48 different sky regions (delineated in HEALPix 
format at A'sido = 2 resolution). 

This resolution is chosen as a compromise between (a) the 
assumption of uniform polarisation fraction in any given sky-region 
and (b) the need to average a large number of data points in order 
to minimise uncertainties. 

The anomalous low-frequency emission correlated with 
thermal-dust is clearly detected in intensity over the entire sky at 
the WMAP K- and Ka- bands. It is also found to be the dominant 
foreground at low Galactic latitude in the pixel analysis, between 
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Figure 10. Histograms of the free-free polarisation fraction at K-band for tliree regions, from 500 realisations at K-band, for the case with no dust polarisation 
(red line) and with 3.5 % dust polarisation (blue line). The input values are shown as vertical dotted lines. The green line and error-bar represents the best-fitting 
value and its Itr uncertainty from the real data. 



b — —40° and b — 10°. The synchrotron spectral-index obtained 
from the K- and Ka-bands is consistent with previous results, al- 
though small differences are found between the polarisation and 
temperature results in the pixel-analysis. On the other hand, the 
values from the all-sky analysis (/? = —3.32 ± 0.12 for intensity 
and /3 — —3.01 ± 0.03 from polarisation) are compatible with an 
average value of /? ~ — 3. 

The polarisation fraction of the synchrotron is constant in fre- 
quency and increases with latitude from 5% near the Galactic plane 
to about 40% in some regions at high latitude; this is consistent 
with a depolarisation effect due to integrating through high-column 
densities in the plane. The average synchrotron fractional polarisa- 
tion at low latitude (|6i < 20°) is 8.6 ± 1.7(stat) ± 0.5(sys)%, 
while at high latitude (16j > 20°) is 19.3 ± 0.8 (stat) ± 0.5 (sys) % 
with a standard deviation of ±11.0%. The robustness of these val- 
ues was verified using monte carlo simulations. 

Previous work has usually assumed that the anomalous dust- 
correlated emission in K-band is unpolarised. We find this as- 
sumption to be reasonable with most of the sky having polarisa- 
tion fractions less than 10 %. However, monte carlo simulations 
revealed that some pixels are likely to be affected by cross-talk, 
at up to the ~ 5 % level. Taking an average of the significant 
(> 2a) regions led to a bias of « 1.5% in the simulations. 
The average polarisation fraction, for the more reliable pixels, was 
3.2 ± 0.9 (stat) ± 1.5 (sys) %. Given the average bias seen in the 
simulations (1.5 %), we do not claim a strong detection, while it is 
clear that the polarisation fraction is less than 5 % at 95 % confi- 
dence. We found several pixels that appeared to have a significant 
detection of dust polarisation at K-band while being robust in the 
simulations. Eight out of 48 regions exceeded the 99th percentile 
of the distribution from simulations with no input foreground po- 
larisation, 6 of which are detected at above 2a and display po- 
larisation fractions between 2.6 % and 7.2 %, except for region 
45 at {l,b) ~ (46°, —66°), which was found to be polarised at 
31.7 ± 11.7 %(2.7o-). 

Free-free emission is found to be unpolarised over the entire 
sky with an upper limit on the fractional polarisation of about 5% 
in nine regions, and about 10% in seven regions. Unlike for the 
dust, no regions were found to have a polarisation fraction that ex- 
ceeded the simulation distributions. Guided by simulations, we find 
an upper limit of 3.4 % (95 % confidence). 
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Table Al. Full fit temperature coefficients for dust in all 48 pixel regions, Table A2. FuU fit coefficients for polarized dust in aU 48 pixel regions, witii 

witli associated 1-a errors (in fiKfiKp^a). associated I-ct errors (in fiKfj,Kpjy„). 
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A T ±1 

U.2 


4 


o 4 ±0.9 


4 4 ±0.9 


1 /I ±0.9 
3.4 -^^--^ 


s ±1 
2.8 


1 q ±1.3 


4 


±0.2 


A ^; ±0.2 
0.0 -^^•■^ 


n o ±0.2 
0.2 -^^-^ 


A 1 ±0.3 

-0.1 -^^-^ 


n T ±0.3 
-U.2 


5 


6 


1 ±0.4 


1 ±0.4 
2.3 -^^-^ 


1 Q ±0.4 

I.y -^^-^ 


1.5 ±0.5 


5 


1 /; ±0.1 
1.0 -^^-^ 


A ±0.1 
0.5 -^^-^ 


A ±0.1 
0.3 -^^-^ 


A 1 ±0.1 

0.1 -^^-^ 


A 1 ±0.1 
U.l -^'^■^ 


g 


0.6 


1 ±1.6 
-1.2 -^^-^ 


1 T ±1.5 
-1.7 -^^-^ 


T ±1.8 


-1.7 ^^-^ 


5 


A < ±0.3 
0.5 -^"-^ 


A /I ±0.3 
0.4 -^"-^ 


A ±0.3 


A ±0.3 


A 1 ±0.4 
0.1 -^"-^ 


7 


A 8 ±1-2 

f .O 


±1.3 


1 c ±1.5 
1.5 -^^-^ 


n a ±1.6 
0.6 -^^-^ 


7 ±1.6 


7 
1 


A ±0.3 

o.y -^"-^ 


A ±0.3 
0.2 -^"^-^ 


A 1 ±0.3 
-0.1 -^"-^ 


A ±0.3 


A ±0.4 
0.2 -^^-^ 


g 


2.3 ±1 


1 ±11 

-i.y 


Q ±1-2 


^ 1 ±1.2 


.3±i 


g 


1 ±0.2 
1 


A ±0.2 
0.2 


A ±0.2 



A 1 ±0.3 
-U.l 


A 1 ±0.3 
0.1 


Q 


4.2 ^l^i-* 


A < ±1.9 

0.5 


A 1 ±2 

-0.1 


f\ <; ±2 
-0.5 


n 9 ±2-4 


Q 


A n ±0.5 
0. / -^"-^ 


A < ±0.5 
0.5 


A -2 ±0.5 

0.3 -^"^-^ 


A 1 ±0.6 
0.1 -^"^ 


A -2 ±0.8 

0.3 


10 


6.5 ii''^ 


T 1 ±1.8 


1 1 ±1.8 
2. 1 


1 /I ±1.9 
1.4 ^^-^ 


1.2 ±1-8 


10 


1 ±0.4 
1.2 -^'^■^ 


A ±0.4 


A ±0.4 


A ±0.4 


A ±0.5 


[ I 


5.4 ±^ 


1 ^ ±1 
1.5 


n o ±1 
0.2 


n 1 ±1.2 


-0.4 ±1 


1 1 


1 ±0.3 
1.2 


A ±0.3 


A 1 ±0.3 
-0.1 


A ±0.3 
0.2 


A ±0.4 
0.2 -^'^■^ 


12 


5.9 ±1 


Q ±1 

2.y 


T ±1 

2 


1 5 ±1.2 
1.5 ^^-^ 


J ±1.3 


12 


1 1 ±0.3 
1.3 -^'^■'^ 


n ^; ±0.3 
0.0 -^^■'^ 


A T ±0.3 

0.2 -^'-'■^ 


A 1 ±0.3 
0.1 -^'^■'^ 


A 1 ±0.4 

0.1 -^^-^ 




c , ±1.1 

J. i 


n ±1-1 
U.2 -^^-^ 


1 1 ±1.2 
-1.1 -L^-^ 


1 ^ ±1.1 


-1.5 


13 


T 1 ±0.4 
2.3 -^''-^ 


n T ±0.4 


A 1 ±0.4 
-0.1 -^"^-^ 


n -3 ±0.5 
-0.3 -^'-'■^ 


n Q ±0.6 
0.3 -^'-'■'J 


14 


y 2 ±0-3 


1 A ±0.3 


2.5 


1 o ±0.4 
1.9 ^^'-^ 


1.4 


14 


1 T ±0.1 


^^ z ±0.1 
0.5 -^'-'■^ 


n a ±0.1 
0.3 -^'^'■^ 


A 1 ±0.1 

0.1 -^''--^ 


A ±0.1 


1 < 

i J 


-i.D 


/I ±0.9 


/I -7 ±1 

-4.7 


A a ±11 
-4.9 -^^-^ 


-4.2 ±1-1 


1 ^ 

1 J 


1 ±0.2 
1.0 -^^-^ 


A A ±0.2 
0.0 -^^-^ 


A ±0.1 
0.2 -^^-^ 


A ±0.2 


A 1 ±0.2 
-0.1 


16 


22 9 ±0.8 


-7 Q ±0.9 

/.y 


A Q ±0.9 

u.y 


< 8 ±0.9 
5. is 


5.1 ±1-2 


16 


A ±0.2 
2.0 


A -7 ±0.2 
0. / 


A a ±0.2 
0.3 


A 1 ±0.2 
0.1 


A ±0.3 
0.2 


17 


<: <: ±0.9 

J.J 


-2 ±1 

2.3 


1 -7 ±1.1 

1.7 -^^-^ 


1 A ±1.2 
1.4 


8 ±1 


17 


1 /: ±0.3 
l.u -^"-^ 


A < ±0.3 
0.5 -^"^-^ 


A 1 ±0.3 
-0.1 -^"-^ 


A ±0.3 


A 1 ±0.4 
0.1 -^'^■^ 


18 


1.2 


1 ±1-4 
-1.2 


1 T ±1.4 

-1. / 


1 « ±1.4 
-1.0 


-1.9 ^i-s 


18 


1 ±0.3 
1 


A ±0.3 



A /I ±0.3 
0.4 


A ±0.4 

u 


A -2 ±0.5 

0.3 


19 


3 =•=1-1 


n ±1-1 

-u.y 


1 < ±1.2 

-1.5 -^^-^ 


1 ±1.4 
-1.8 -^^-^ 


-1.7 ±^■2 


19 


1 1 ±0.3 
1.1 -^^-^ 


A ±0.3 
0.2 -^"-^ 


A ±0.3 
0.2 


A ±0.4 


A 1 ±0.4 
-0.1 -^"-^ 


20 


9 ] ±0.3 


-5 ±0.3 

3.y 


^ ±0.3 
2.0 


1 £ ±0.3 
1.0 


1.1 


20 


1 -J ±0.1 
1.3 


A 'X ±0.1 
0.3 


A 1 ±0.1 
0.1 


A ±0.1 




A ±0.1 




21 


-J ±0.7 


a 1 ±0.7 
3.1 


T T ±0.8 

2.2 


1 1 ±0.8 
1.1 ^''-^ 


1.1 ±1 


21 


A 1 ±0.3 
-0.1 -^'^■'^ 


A 1 ±0.3 
0.1 -^'^■'^ 


A ±0.3 


A ±0.4 


A /I ±0.5 
0.4 -^'^■'^ 


22 


8.3 ±^'2 


±1.2 


c 1 ±1.3 
5.1 


/ T ±1.2 


3.5 


22 


A ^ ±0.3 
0.5 


A ±0.3 



A 1 ±0.2 
-0.1 


A ±0.3 



A 1 ±0.4 
-0.1 


23 


9.9 ±1'' 


1 ±1.7 
5.3 


1 Q ±1.8 
3.0 -^^■'^ 


T 1 ±1.9 

3.1 


2.4 ±2.3 


23 


1 1 ±0.5 


n a ±0.5 
0.3 


n 1 ±0.5 
0.1 -^'^'■'^ 


n a ±0.6 
0.3 -^''-'^ 


n 1 ±0.8 


24 


3.8 


1 « ±2 


1 « ±2.2 
1.0 


1 ^ ±1.9 
1.3 


0.3 ±2-5 


24 


1 -7 ±0.6 
1. / 


A ±0.6 
0.2 


A /I ±0.6 
0.4 


A A ±0.8 

-O.D 


A « ±0.9 
0.0 


25 


.3 7 ±0.9 


-7 ^; ±0.9 
-7.D 


Q /; ±0.9 
-0.0 -^^-^ 


« Q ±1.1 
-5.9 


-7.8 ±^■2 


25 


1 c ±0.3 
1.5 -^^'^ 


A 1 ±0.3 


A ±0.3 
0.2 -^^■'^ 


A 1 ±0.3 
0.3 -^^-^ 


A ±0.4 


26 


2.6 ^i-''' 


1 1 ±1.9 
3.1 -^^-^ 


1 /I ±1.9 
3.4 -^^-^ 


-J /I ±2 
3.4 


3.4 =■=2.4 


26 


A A ±0.4 

0.0 -^^-^ 


A -J ±0.4 
0.3 -^'^■^ 


A 1 ±0.4 
U.l -^"^-^ 


A 1 ±0.5 
0.1 -^"-^ 


A 1 ±0.6 
0.1 


27 


20 4 ±0.8 


A 1 ±0.9 
0.1 -^'^■^ 


c ±0.9 
5 


■3 ±0.9 
3.8 


2.9 ±1 


27 


^ ±0.2 


A ±0.2 
0.2 -^"■'^ 


A 1 ±0.2 

-0.1 -^^-^ 


A 1 ±0.3 
-0.1 -^"^-^ 


A ±0.3 


28 


11.3 ±2.2 


±2.4 

D.y -^-^-^ 


c 1 ±2.6 
5.1 -^^■'^ 


A 1 ±2.6 
4.1 ^^-"^ 


3 ±2-6 


28 


n ±1-1 
U.5 -^'-■■^ 


A ±11 
-0.2 -^'-■■^ 


1 ±1 
1 


A ±1.3 
0.2 -^^■'^ 


J ±1.5 


29 


y 2 ±0-6 


1 -7 ±0.6 
1.7 -^'^■'J 


A 1 ±0.7 
U.3 


n 7 ±0.7 
-0.7 


-0.6 ±° '7 


29 


1 /I ±0.2 
1 .4 -^'-'■^ 


n c ±0.2 
0.5 -^'-'■^ 


A o ±0.2 
0.2 -^'^■^ 


n ±0.2 


A ±0.2 


30 


4.2 


n ±0.6 

U.y -^^'-^ 


A T ±0.7 

U.2 


n 1 ±0.6 


.0.3 ±0.7 


30 


n 1 ±0.1 
U.3 -^'■'■-^ 


A 1 ±0.2 
0.1 -^^-^ 


A ±0.2 


A ±0.2 


A ±0.2 


31 


-0.6 


1 ^ ±1-2 

-3.5 


a ±1.3 


4 ±1.2 


-3.9 ±1-3 


31 


1 3 ±0.2 
1.3 -^'-'■^ 


n 1 ±0.2 
0.3 -^"-^ 


A 1 ±0.2 
0.1 -^"-^ 


A 1 ±0.3 
-0.1 


A 1 ±0.3 
0.1 -^^■'^ 


32 


7.2 ±1 


a ±1-1 

3.y -^^-^ 


-3 ±1.1 


Z.J ' 


1.7 ±1-3 


32 


T A ±0.2 
2.D -^'■■'■^ 


1 ±0.2 


A A ±0.2 

0.4 -^'-'■^ 


A 1 ±0.3 
-0.1 


A ±0.3 
-0.2 -^^-^ 


33 


4.6 ±''■3 


n /I ±0.9 
U.4 -^'-'■■^ 


n /I ±0.9 


1 ±0.9 


-1.3 ±1 


33 


o c ±0.3 
2.5 


n ±0.3 
0.0 -^^■'^ 


A ±0.3 
0.2 -^'^■'^ 


n o ±0.4 
0.2 -^"^-^ 


A /I ±0.5 
-0.4 


34 


14.6 ±1-3 


in T ±1.3 
1U.2 -^^-^ 


Q -7 ±1.4 


7 1 ±1.6 

/.I 


6 ^i-"" 


34 


]^ ±0.3 


A 1 ±0.3 
-0.1 -^^"^ 


A 1 ±0.3 
0.1 -^'^■'^ 


A 1 ±0.4 

-0.1 -^^-^ 


A 1 ±0.5 
-0.1 


35 


^i i^ 


-7 ±1.7 
-2.7 


1 1 ±1.6 
-3.1 -^^-^ 


1 ±2 
-3 


-2.8 * 


35 


A A ±0.3 

O.u -^^-^ 


A 'J ±0.3 
0.3 -^"-^ 


A ±0.3 


A ±0.4 
-0.2 -^"-^ 


A < ±0.5 
0.5 -^"-^ 


36 


4 4 ±0.5 


A 1 ±0.5 


1 ±0.6 
-1.2 -^"^-^ 


1 Q ±0.7 

-i.y 


-1.8 


36 


1 ±0.1 
1.2 -^"-^ 


A /I ±0.1 
0.4 -^"-'■^ 


A ±0.1 

0.2 -^^-^ 


A 1 ±0.1 
0.1 -^"-^ 


A 1 ±0.2 
-0.1 -^"-^ 


37 


7.1 ±1-3 


1 1 ±1.4 

3.1 -^-^-^ 


1 ±1.3 
I.y -^^■'^ 


2 ±1.3 


1.1 


37 


1 < ±0.4 
1.5 -^''-^ 


A 1 ±0.4 
0.3 -^"-^ 


A a ±0.4 
0.3 -^'-'■^ 


A 1 ±0.4 
-0.1 -^"-^ 


A 1 ±0.5 
-0.3 -^"-^ 


38 


-1.2 ii-'- 


-3 A ±1 1 


a /I ±1 2 
-3.4 =^-L-^ 


Q c ±1 1 


-3 ±1-1 


38 


/I ±0 3 
0.5 ^^'-^ 


A A ±0 3 
0.0 -^'-'■^ 


A ±0 3 


A 1 ±0 4 

-0.1 =^'~'-^ 


A ±0 4 


39 


5.8 ±2 


n ±2.1 


A -7 ±2.2 

u. / -^^-^ 


-0.4 ±2-3 


0.3 ±2.8 


39 


A -7 ±0.6 


A 1 ±0.7 
0.3 


A 1 ±0.6 
0.1 -^'^■'^ 


A 1 ±0.8 
0.1 


A 1 ±0.9 
-0.1 -^"-^ 


40 


-6.6 ±1-2 


1 n -7 ±1.3 


1 1 1 ±1.3 
-1 1.3 -^-^-^ 


11 1 ±1.3 
-1 1.3 ^^-^ 


-9.5 ±1-5 


40 


A ±0.3 
0.5 -^^-"^ 


A a ±0.3 
0.3 -^^■'^ 


A 1 ±0.3 
0.1 -^'^■'^ 


A ±0.3 


A ±0.4 


41 


-3 ±1'5 


-5 =^1-^ 


-5.8 ±1-7 


-6.1 ±1-6 


4 ±1.8 


41 


0.3 ±"-4 


0.1 ±"-4 


-0.1 ±0-4 


±0.5 


-0.2 ±0.6 


42 


2.4 ±2-2 


-0.3 ±2.3 


-1 ±2.2 


-1.4 ±2-7 


-1.5 =^2.9 


42 


1.8 ±^-5 


4 ±0.5 


0.3 ±0.4 


0.1 ±0-6 


0.1 ±0-7 


43 


-0.1 ±1-3 


-2.3 ±1-4 


-2.9 il-""' 


-3.2 ±1-4 


-2.6 ^^-^ 


43 


I ±0.2 


0.5 ±"-2 


0.1 ±0-2 


-0.1 ±0-2 


-0.1 ±0-3 


44 


7.7 ±1 


2.9 ±1 


1.3 ±1 


0.6 ±11 


0.4 ±1-1 


44 


3 ±0-2 


0.9 ±0-2 


0.6 ±0-2 


-0.1 ±0-3 


0±0.3 


45 


6.1 ±1-8 


3.1 ±1-5 


2.6 


2.2 ±1-8 


2 ^^-^ 


45 


1.5 ±0-4 


0.9 ±0-4 


0±0.4 


0±0.5 


0.1 ±0-6 


46 


12.7 ±1-8 


8.9 ±1-8 


8 ±2.1 


7.7 ±2 


6.1 ±1-9 


46 


0.1 ±0-5 


0±0.5 


0.2 ±0-5 


-0.1 ±0-6 


-0.6 io '- 


47 


5.6 ±2 


0.1 ±2.1 


-0.9 ±2-4 


-2 ±2.2 


-1.8 ±2.7 


47 


0.8 ±0-4 


0.4 ±0-4 


0.7 ±0-4 


0.2 ±0-5 


-0.1 ±0-6 


48 


.1.4 ±1.8 


_4 ±2 


-4.4 ±1-9 


-4.8 ±2.2 


-3.9 ±2.4 


48 


2 ±0.4 


0.7 ±0-4 


0.4 ±0-3 


-0.2 ±0-4 


0±0.6 



Table A3. Full fit temperature coefficients for Synchrotron in all 48 pixel Table A4. Full fit coefficients for polarized synchrotron in aU 48 pixel re- 

regions, with associated \-a errors (in M^^408MHz-'- gions, with associated \-a errors (in M^^408mhz)- 
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Pixel K-band Ka-band Q-band V-band \\ -band Pixel K-band Ka-band Q-band V-band VV-band 

Free-Free Temperature C-C coefficients Free-Free Polarization C-C coefficients 



1 
1 


<Q Q ±23.9 


A A "7 ±25.1 
-44. / 




07 ±29.6 

-JZ ^ " 


77 1 ±33.4 


1 
1 


7 ±5.8 

-J -^^ ^ 


n e ±5.9 

u.o 


<; 4 ±5.6 
-J. 4 


n 7 ±6.8 
-U. / 


1 <; ±8.3 
1 .J 


z 


£7 o ±21.7 
O/.Z 


77 ^ ±22.8 
/ / .O 


ori Q ±26.1 

ou.y 


7Q ±27.6 

/y 


fi7 ±23.5 
oz.o 




z 


7 ±5.7 
-Z.O 


n ^ ±5.7 

U.J 


1 7 ±5.4 


1 ft ±fi-*' 

1.0 


7 Q ±7.9 

/.y 


J 


S ±11-3 
-44. J 


/in A ±11.6 

-4U.D 


/IT 1 ±11.5 


/II 1 ±12.3 
-41.1 


a ±11.8 
-J J. J 


; -3 
J 


in ft ±3.5 

-lU.D 


n 7 ±3.7 

-U.J 


n 1 ±3.6 
-u. 1 


1 ±4.5 
Z.J 


1 ±4.9 
1.0 


A 

4 


00 1 ±3.5 
ZZ.l 


10 ±3.6 
IZ.o 


y.o 


5 7 ±4.6 
J. / 


7 c ±4.8 

Z.J 


4 


1 S ±1.2 
-l.J 


n 4 ±1.3 
-U.4 


n A ±1.2 

-u.4 


n ±1-5 
u 


1 Q ±1.8 

-I.y 


c 

J 


lO 1 ±10.1 
IZ. 1 


O A ±11.1 
(5.4 


in ±10.5 
lU 


/l Q ±12.5 
4.y 


7 A ±11.3 
J. 4 


c 

J 


1 1 7 ±2.8 
-1 l.Z 


n 7 ±2.8 
-u. / 


7 7 ±2.8 
-Z.J 


1 A ±3.4 
-1.4 


4 s ±4.1 
4. J 


D 


o 1 ±8.5 
-5.1 


in /I ±9.2 
-1U.4 


im ±9.7 


1 1 ±9.7 
-1 1.0 


Q 7 ±11.6 

-y.z 





^ ±2 

J.y 


no ±2 
u.z 


n 7 ±1.9 

-U.J 


n 1 ±2.3 
u.l 


n 8 ±2.7 
u.o 


1 


1 /.J 


in A ±8.1 

1U.4 


7 1 ±8.4 
/.I 


c 7 ±8.7 
J. / 


4 e. ±10.3 
4.0 


/ 


4 7 ±1.6 
-4.Z 


1 7 ±1.7 
-1. / 


n 7 ±1.5 

-U.J 


n ±2 
U.Z 


n 7 ±2.4 
-u. / 


Q 

o 


in 1 ±13.1 


in 7 ±14.2 


in ±16.1 

-lU 


Q ±14.2 

-y.y 


<: /; ±14.1 
-J.O 


Q 


9 ±3.6 


1 4 ±3.6 
1.4 


n 4 ±3.4 

-U.4 


n c ±4.2 
-U.J 


n <; ±5 

-U.J 


Q 


o o ±5.4 
-J.J 


7 A ±5.8 
- /.4 


8 1 ±6.7 
-cS.l 


Q ±6.1 
-y.y 


±6.7 
-0 


n 

y 


1 f. ±1.6 
1.0 


n 4 ±1.7 

u.4 


n 4 ±1.6 

u.4 


n 7 ±1.9 
U. / 


1 7 ±2.3 
l.J 




10 /; ±4.8 


1/; f, ±5.1 
10. 


ic -1 ±5.3 
1 J.J 


1 1 7 ±5.9 
11./ 


7 ±5.5 
0. / 




n 7 ±1.2 
-U.J 


1 7 ±1.2 
l.J 


n 1 ±1.2 
u. 1 


n A ±1-4 

-U.4 


1 7 ±1.7 
-l.J 


1 1 


10 ±3.3 

iy.z 


1^; o ±3.4 
10. ?5 


IS ±3.9 

1 J.O 


1-! 7 ±3.9 
1 J. / 


17 f. ±3.4 
iZ.O 


1 1 


n <: ±1.1 
U.J 


n 7 ±1.1 
u.z 


1 ±1.1 
1 


n ±1-4 
u.z 


n 7 ±1.6 
-U.J 


1 


D.4 


n /I ±5 
U.4 


o T. ±5.4 
-Z.J 


£ 7 ±4.9 

-O.Z 


c 7 ±5.5 
-J. / 


1 

iZ 


n ±1.5 

u.o 


n ±1.5 
u 


n 7 ±1.6 
u. / 


n ft ±1-9 

-u.o 


1 Q ±2.3 

-I.y 


1 1 


00 A ±2.1 

zz.o 


1 A Q ±2.2 

14. J 


11 Q ±2.5 
i 1 .5 


7 ±2.2 
/.o 


£ 7 ±2.5 

D.J 


1 1 
1 J 


1 1 ±0.8 
-i.i 


n ±0.8 
u 


n 5 ±0.8 

-U.J 


n 7 ±0.9 

-U. / 


n 7 ±11 
-U.J 


1 /I 
14 


n ±3.7 


-J.Z 


A ±4.3 
-4 


<; 7 ±4.2 
-J.Z 


<; A ±4 
-J. 4 


1 A 
14 


A 7 ±0.8 
-4.Z 


n 9. ±0.8 
-u.o 


1 7 ±0.8 

-l.Z 


n 1 ±0.9 

-U. 1 


-i.i 


1 c 
ID 


ic ±1-4 
Ij.Z 


in 1 ±1-5 

lU.j 


S A ±1.5 


£ 7 ±1.5 
0. / 


c 7 ±1.8 
J. / 


1 c 
Ij 


1 S ±0.4 
-l.J 


n 4 ±0.3 

-u.4 


n 7 ±0.3 
-u.z 


n ±0.4 
u 


n ±0.5 
u 


1 A 


7 ±1.7 

/.z 


o 1 ±1.7 
J.l 


1. / 


n 1 ±1.8 

-U.J 


n s ±2 
-U.o 


10 


n ^ ±0.3 

U.J 


n 1 ±0.3 

U. 1 


n4 ±0.4 
u.4 


n ±0.4 
U 
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Table AS. FuU fit temperature coefficients for free-free in all 48 pixel re- Table A6. Full fit coefficients for polarized free-free in all 48 pixel regions, 

gions, with associated l-cr errors (in ij,KR~^). wiHi associated 1-a errors (in fiKR~^). 
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